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Abstract 

A relativistic mean-field model is used to study the ground-state properties of neutron-rich nuclei. 
Nonlinear isoscalar-isovector terms, unconstrained by present day phenomenology, are added to the 
model Lagrangian in order to modify the poorly known density dependence of the symmetry energy. 
These new terms soften the symmetry energy and reshape the theoretical neutron drip line without 
compromising the agreement with existing ground-state information. A strong correlation between 
the neutron radius of 208 p]^ ^^^^ ^.j^g binding energy of valence orbitals is found: the smaller the 
neutron radius of ^''^Pb, the weaker the binding energy of the last occupied neutron orbital. Thus, 
models with the softest symmetry energy are the first ones to drip neutrons. Further, in anticipation 
of the upcoming one-percent measurement of the neutron radius of ^'^^Pb at the Thomas Jefferson 
Laboratory, a close relationship between the neutron radius of ^'^'^Pb and neutron radii of elements 
of relevance to atomic parity- violating experiments is established. 
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I. INTRODUCTION 



Core-collapse supernovae and their remnant neutron stars are an important source of 
information pertaining to both superdense matter and nuclei far from the valley of stability. 
At present, however, important open questions remain in the understanding of the structure 
and dynamics of neutron stars, such as their size, composition, and cooling mechanism. 
Increased knowledge of these issues may in turn lead to a better understanding of stellar 
burning and heavy-element nucleosynthesis. 

The structure of spherical neutron stars in hydrostatic equilibrium, the so-called 
Schwarszchild stars, is solely determined by the equation of state of neutron-rich matter 
in beta equilibrium. Such an equation of state is well represented as the sum of two distinct 
components: i) a symmetric {N = Z) component that is well constrained at saturation den- 
sity plus a) a. symmetry energy that accounts for any possible neutron-proton imbalance. 
Some of the important neutron-star properties mentioned above depend critically on the 
poorly known density dependence of the symmetry energy. Indeed, as a result of the satu- 
ration of symmetric nuclear matter, most of the pressure supporting the star up to about 
two times saturation density is provided by the symmetry energy. 

In a set of recent papers, the sensitivity of the symmetry energy to changes in the model 
Lagrangian was investigated [J, 0, 0, ■ The model Lagrangian was modified by introducing 
nonlinear couplings between the isoscalar and the isovector mesons. These new terms enable 
one to modify the neutron skin of heavy nuclei without changing ground-state properties 
that are well constrained experimentally. It should be noted that precise information on 
the neutron radius of a heavy nucleus — 208p]^ — might soon become available via a parity- 
violating electron scattering experiment at the Jefferson Laboratory that promises a 1% 
accuracy |^, |^ . While such a measurement will provide the most accurate determination 
of the neutron density of a nucleus to date, it will also impact strongly on astrophysical 
observables. Indeed, a model-independent (or data-to-data) relation between the neutron 
skin of ^'^spb and the crust of a neutron star was recently established This strong 
correlation emerges as a result of the similar composition of the neutron skin of a heavy 
nucleus and the crust of a neutron star, namely, neutron-rich matter at similar densities. 
Further, the measurement of the neutron skin of ^ospb constrains the cooling mechanism in 
neutron stars. A small neutron radius {Rn) suggests a soft symmetry energy which favors 
a small proton fraction in dense matter. In turn, a low proton concentration rules out the 
enhanced cooling of neutron stars via the direct URCA process Q. If URCA cooling is 
indeed ruled out, then observations of enhanced cooling may provide strong evidence in 
support of exotic states of matter at the core of neutron stars (see, for example, Ref. 0]). 

The density dependence of the symmetry energy should also play an important role on 
the properties of nuclei far from the valley of stability. Indeed, some theoretical calculations 
predict the emergence of new magic numbers as a result of a reduced spin-orbit splitting 
originating from a diffuse neutron skin |^. No place is this more evident than in the 
neutron-rich oxygen isotopes. The nucleus of ^^O, which in the traditional scheme has filled 
1^5/2 2s^l2 orbitals, appears to be the heaviest member of the isotopic chain that remains 
stable against particle emission. Yet, most theoretical calculations predict the existence of 
the "doubly magic" nucleus ^^O. This issue continues to be revisited in light of a remarkable 
experiment that shows that the mere addition of a single proton stabilizes the fluorine chain 
up to ^^F; this is six more neutrons than in ^^O [lo| . 

In this — our first — study of neutron-rich nuclei, we focus on the effect of the new isoscalar- 
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isovector terms on ground-state properties. This study will be conducted within the frame- 
work of relativistic mean-field models that reproduce a large body of ground-state observ- 
ables for nuclei at, or near, the valley of stability. Yet these observables are mostly insensitive 
to the density dependence of the symmetry energy. It is the aim of the present study to 
constrain the symmetry energy through a systematic study of ground-state properties of 
neutron-rich nuclei. While we recognize that additional physics, such as pairing, may be 
needed for more realistic studies, we limit ourselves to study the impact of the bulk sym- 
metry energy on the ground-state properties of neutron-rich nuclei. Further, we focus on 
self-consistent mean-field models in the hope of extending earlier studies of the linear re- 
sponse jn| to incorporate the new isoscalar-isovector terms. 

The manuscript has been organized as follows. In Sec. |n]we review briefly the relativistic 
formalism paying special attention to the role of the new isoscalar-isovector couplings. We 
illustrate how these new couplings modify the density dependence of the symmetry energy 
while leaving unchanged all properties of symmetric nuclear matter. In Sec. Illll we present 
the results of the calculations for various ground-state observables for a variety of nuclei. A 
summary and conclusions are presented in Sec. IIVI 



II. FORMALISM 

In this section we describe in some detail the model Lagrangian employed in this work 
and the mean-field approximation used to compute ground-state properties for a variety of 
neutron-rich nuclei. Although most of the derivations are standard, a review is included 
here to illustrate the role of the mixed isoscalar-isovector meson terms first introduced in 
Ref. These terms, which supplement the phenomenologically successful Lagrangian of 
Ref. modify the density dependence of the symmetry energy, thereby reshaping the 
nuclear landscape away from the valley of stability. 



A. Effective Lagrangian 



The Lagrangian of Ref. jlj] includes an isodoublet nucleon field {ijj) interacting via the 
exchange of one scalar (0 for the sigma) and three vector {V^ for the omega, b'^ for the rho, 
and A^^ for the photon) fields. That is. 
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where the various field tensors have been defined as follows: 
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In addition to meson-mediated interactions, the Lagrangian is supplemented with nonlinear 
meson interactions that serve to simulate the complicated dynamics that lie beyond the 
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realm of the low-energy effective theory. Indeed, by fitting the constants that parametrize 
these meson terms to the bulk properties of nuclei, rather than to two-nucleon data, the 
complicated dynamics originating from nucleon exchange, short-range effects, and many- 
body correlations get implicitly encoded in a small number of empirical constants. For the 
purpose of the present discussion, the following local meson terms are sufficient: 

?7eff(0, b'^) = + - ^{W.W^'f - (As$2 + KW.W^') (B^ ■ B'^) , (3) 

where the following definitions have been introduced: ^ = gs(p, W^ = g^Vfj,, and B^ = 5fpb^. 
The inclusion of scalar-meson interactions {k and A) is dictated by the empirical value of 
the compression modulus of symmetric nuclear matter at saturation density {K = 200 — 
300 MeV). In contrast, quartic vector self-interactions {() affect primarily the high-density 
component of the equation of state; their impact at densities below normal nuclear-matter 
saturation density is yet to be determined 0|. Finally, the nonlinear mixed isoscalar- 
isovector couplings (As and Av) are powerful terms because they can be used to modify the 
density dependence of the symmetry energy ^ while not influencing ground state properties. 
While power counting suggests that other local meson terms (such as mixed scalar-vector 
cubic terms and quartic rho-meson self-interactions) may be equally important jl^], their 
phenomenological impact has been documented to be small 0, E3) they will be not be 
considered in this study. 



B. Mean-field Equations 

The field equations resulting from the above Lagrangian may be solved exactly in the 
mean-field limi t by replacing all meson-field operators by their expectation values, which are 
classical fields |l3|. For a static, spherically symmetric system this implies (using |x|=r): 



0(x) ^ (0(x)) = Mr) , (4a) 

Vf^ix) ^ {V^{x)) = ^?^Vo(r) , (4b) 

Ki^) ^ m^)) = g'Xsboir) , (4c) 

A>^{x) ^ (A^(x)) = g^'Aoir) . (4d) 

Similarly, the various baryon sources to which these mesons couple must be replaced by 
their (normal-ordered) expectation values in the mean-field ground state. That is, 

ip{x)l^{x) {■.^{x)li){x):) = ps{r) , (5a) 

^(a;)7'^V(a;) (:^(x)7'^^/'(x) :) = g^^ p.{r) , (5b) 

^(a;)7^r,V'(a:) (:^(x)7'^r,V'(x) :) = (?'^°W3(r) , (5c) 

^(x)7%^(x) ^ (:V^(x)7%^(x):) = (7^'Vp(0 • (5d) 

Note that the proton isospin projection operator has been defined as Tp = (l + r3)/2. The 



baryon sources generate (classical) meson fields that satisfy coupled, nonlinear Klein-Gordon 
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equations of the following form: 
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The photon field couples only to the (point) proton density and its solution is thus reduced 
to quadratures 

~ pr POO 

dxx^Pp{x)+ / dxxpp{x) . (7) 
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The eigenstates of the Dirac equation, for the spherically symmetric mean-field ground 
state assumed here, may be classified according to a generalized angular momentum x. 
Thus, the single-particle solutions of the Dirac equation may be written as 
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where C,t denotes a two-component (Pauli) spinor in isospin space (with t = il=l/2 for protons 
and neutrons, respectively), n and m are the principal and magnetic quantum numbers, 
respectively, and the spin-spherical harmonics are defined as 



y>cm{y^) = {yi\l\jm) ; j = |x| 



X , if X > 0; 

-1 — X , if X < 0. 



(9) 



Note that the phase convention adopted in Eq. (jHJ [i.e., the relative factor of i) is such that 
real bound-state wave functions {g and /) are generated if the mean-field potentials are also 
real. Further, the following spinor normalization has been adopted: 
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where {a} = {a = nxt; m} denotes the collection of all quantum numbers required to de- 
scribe the single-particle Dirac spinor. It then follows that the coupled differential equations 
satisfied by the radial components of the Dirac spinor are given by 
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where the upper and lower numbers in these equations correspond to protons and neutrons, 
respectively. Having determined all occupied single-particle states, the various ground-state 
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densities, which act as sources for the meson fields in the Klein-Gordon equations [see 
Eq. ©], may now be computed. They are given by 

ps{r) = Ps,p(r) + ps,n(r) , (12a) 

Pvir) = pv,p{'^) + PvA'^) ' (12b) 

P3(r) = Pv,p(r) -pv,n(r) , (12c) 

Pp(r) = Pv,p(r) , (12d) 

where scalar and vector densities have been defined as 

Ground-state properties of the system are described by a solution of the coupled, ordinary 
differential equations for the classical meson fields [Eq. ^] and for the Dirac single-particle 
states [Eq. PHI ]. The solution must be self-consistent, that is, the meson fields generating 
the Dirac mean-field potentials must satisfy Klein-Gordon equations having ground-state 
densities constructed from the same single-particle states as their sources. Thus, an it- 
erative procedure must be implemented. The self-consistent procedure starts with initial 
Woods-Saxon shaped meson fields of reasonable strength and range to generate, via a con- 
ventional Runge-Kutta algorithm, bound-state energies and corresponding wave functions 
for all occupied single-particle states. At this point, scalar and vector densities (ps,t and 
Pv,t) are computed. Using these as sources for the meson-field equations, new meson fields 
are generated by using Green's function techniques. The newly generated meson fields will 
differ, in general, from the initial Woods-Saxon guess. Thus, this iterative procedure must 
continue until self-consistency (convergence) is achieved. 



III. RESULTS 

The symmetry energy of infinite nuclear matter impacts on the dynamics of neutron-rich 
nuclei. After all, the symmetry energy describes how the energy of nuclear matter increases 
as the system departs from equal numbers of neutrons and protons. To investigate the 
structure of neutron-rich nuclei we use a variety of effective field theory (EFT) models that 
differ in their prediction for the density dependence of the symmetry ene rgy . 

Three models will be considered in this text: the very successful NL3 0,0] along with 
the newer S271 and Z271 [1] parameter sets. The S271 and Z271 models are alike the NL3 
in that they are constrained to the following properties of symmetric nuclear matter: i) 
nuclear saturation at a Fermi momentum of kp = 1.30 fm~^, ii) a binding energy per nucleon 
of 16.24 MeV, and in) a compression modulus of = 271 MeV. The first table (Table H} 
lists the various parameter sets that are needed to reproduce these properties of symmetric 
nuclear matter at the mean-field level. 

In Fig. n the equation of state for symmetric nuclear matter is displayed for the three 
models discussed in the text. That all models are identical at (and near) saturation den- 
sity follows from the fitting procedure that has produced the above mentioned constraints. 
While significant discrepancies among the models emerge at high density {kp^l.bO fm~^), 
presumably due to differences in the input values for the effective nucleon mass M* and 
(, these discrepancies disappear at the densities (fcp ^ 1-30 fm~^) relevant to the physics 
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TABLE I: Model parameters used in the calculations. The parameter k and the scalar mass are 
given in MeV. The nucleon, rho, and omega masses are kept fixed at M = 939, rrip = 763, and m^j = 
783 MeV, respectively — except in the case of the NL3 model where it is fixed at m^ = 7S2.f) MeV. 



Model 


rus 
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9l 
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K 


A 
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NL3 


508.194 


104.3871 


165.5854 


79.6000 


3.8599 


-0.01591 


0.00 


S271 


505.000 


81.1071 


116.7655 


85.4357 


6.6834 


-0.01580 


0.00 


Z271 


465.000 


49.4401 


70.6689 


90.2110 


6.1696 


+0.15634 


0.06 



of finite nuclei (see inset in the figure). Thus, the equation of state for symmetric nuclear 
matter is model independent in this range and the only bulk property of infinite nuclear 
matter that can lead to a model dependence in ground-state observables of finite nuclei is 
the symmetry energy. 




0.5 1 1.5 2 2.5 



kp(fm"') 

FIG. 1: Equation of state of symmetric nuclear matter for the three models considered in the text. 

Unfortunately, the density dependence of the symmetry energy is poorly known. Indeed, 
even the symmetry energy at saturation density is not well constrained experimentally. It is 
some average between the symmetry energy at saturation density and the surface symmetry 
energy that is constrained by the binding energy of nuclei. As a consequence, we adjust the 
value of the NNp coupling constant to reproduce a symmetry energy of asym = 25.67 MeV 
at a Fermi momentum of /ci? = 1.15 fm~^ (or p = 0.10 fm~^) (For a recent discussion 

on the surface symmetry energy see Ref. jig.) 

In Tabl e Ull predictions for the binding energy per nucleon [17], root-mean-square charge 
radius (3j llSj - and neutron radius of ^ospb are displayed for the three models considered in 
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TABLE II: Symmetry energy at saturation density, binding energy per nucleon, root-mean-square 
charge radius (point-proton radius in parenthesis) and neutron radius of ^ospj-, £qj, ^^le various 
models considered in the text (with Av=0). 



Model al%, (MeV) B/A (MeV) RchjRp) (fm) Rg (fm) 

NL3 37.285 7.854 5.509 (5.460) 5.740 

S271 36.637 7.939 5.509 (5.460) 5.714 

Z271 36.298 7.775 5.508 (5.459) 5.700 

Exp ~ 37 7.885 5.504 unknown 



the text. The predictions are within one percent of the experimental values, except for the 
neutron radius which is poorly known. The nonlinear coupling (Av) between the isoscalar 
and the isovector mesons [see Eq. Q] enables one to modify the density dependence of the 
symmetry energy, and thus the neutron radius of ^°^Pb, while leaving well-known ground- 
state properties intact. This suggests that existing ground-state information, such as charge 
densities and binding energies, do not determine the neutron radius uniquely. Thus, a 
new measurement, such as the neutron radius in ^osp^^ jg needed to provide important 
constraints on the density dependence of the symmetry energy. It is therefore the aim of 
this paper to correlate ground-state properties of neutron-rich nuclei to the neutron radius 

of 208pb. 
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FIG. 2: The symmetry energy in the NL3 model. 

The model dependence of the symmetry energy is generated by tuning the isoscalar- 
isovector term. (For simplicity, we set the isoscalar- isovector coupling Ag to zero henceforth.) 
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Note that changing Av has no effect on the bulk properties of symmetric nuclear matter as 
the expectation value of the rho-meson field is identically zero in the N = Z limit. It does 
affect the value of the symmetry energy so a mild adjustment of the NNp coupling constant 
is required to keep the symmetry energy fixed at agym = 25.67 MeV; this ensures that the 
binding energy per nucleon in ^osp]^ remains fixed at 7.87 MeV jl7l |. 

In Fig. 121 the symmetry energy of infinite nuclear matter is plotted using the NL3 param- 
eter set for different values of the isoscalar-isovector term Ay (the original NL3 parametriza- 
tion has Av = ^3])- Increasing A^ softens the symmetry energy thereby reducing the 
internal pressure of the system. As a result, the original NL3 model predicts both the 
largest neutron radius in ^osp]^ ^^^^ ^j^g largest neutron-star radii (for a given mass) 0, 
The inset in the figure shows the symmetry energy at the densities relevant to the dynamics 
of neutron-rich nuclei. Note that because all parameterizations are constrained to have the 
same symmetry energy at kp = 1.15 fm~^, models with a softer symmetry energy have a 
larger symmetry energy at low densities. In these softer models there is a higher price to 
pay for departing from equal numbers of protons and neutrons. This suggests a definite 
correlation: models with the smallest neutron skin (Rn — Rp) in ^osp]^ should be the first 
ones to drip neutrons. 




0.08 



0.06 



0.04 



0.02 



FIG. 3: Neutron and proton densities for ^*^Ca in the NL3 model for various values of the isoscalar- 
isovector term Av The inset shows the surface-peaked nature of the (square of the) upper compo- 
nent of the least bound 1/^/^-neutron. 



Further evidence for this behavior is provided in Fig.Elwhere proton and neutron densities 
for the neutron-rich nucleus ^°Ca are displayed. The nonlinear isoscalar-isovector term Ay 
is used to change the density dependence of the symmetry energy which in turn modifies 
the neutron radius of N ^ Z nuclei. Indeed, the predicted values for the neutron skin of 
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60Ca are i?„-i?p = 0.608, 0.567, 0.523 fm for A^ = 0, 0.015, 0.030, respectively. Note that the 
proton radius remains fixed at Rp = 3.562 ± 0.013 fm. The inset shows the square of the 
upper component of the least bound (1/^'^^) neutron in ^°Ca (for Av = 0) to illustrate how it 
is most sensitive to the symmetry energy in the low-density surface region; this is the region 
(according to the inset in Fig. El) where one expects the model with the largest neutron skin 
to give the strongest binding energy for the If^^'^ neutron. In accordance with this statement 
we obtain a binding energy for the If^^'^ neutron in ^°Ca of 4.907,4.705,4.443 MeV for a 
corresponding neutron skin of Rn — Rp = 0.608, 0.567, 0.523 fm, respectively. 
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FIG. 4: The least bound neutron for various neutron-rich nuclei as a function of the neutron skin 
in ^*^^Pb for the three models considered in the text. A "binding energy" of MeV indicates that 
the single-particle orbital is unbound. 



This correlation is displayed in graphical form in the upper-left-hand panel of Fig. ^ 
where the binding energy of the least bound neutron in ^°Ca is plotted for the NL3 model 
(solid line) as a function of the neutron skin in ^°^Pb. As evinced by the other panels 
in the figure, this correlation holds throughout the periodic table. Moreover, it is model 
independent. Indeed, similar plots have been added for the S271 (dashed line) and Z271 
(dot-dashed line) parameterizations. In all cases the binding energy of the least bound 
neutron increases with increasing neutron skin. Note, however, that the shell ordering is 
not constrained among the models. For ^°Ca the NL3 model predicts the If^^"^ orbital to 
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be the least bound, whereas in the S271 and Z271 models it is the 2p^/^ orbital which is the 
least bound. Moreover, this correlation, namely, that an increase in the isoscalar-isovector 
coupling is responsible for a decrease in both the neutron skin of ^^'^Pb and in the binding 
energy of the least bound neutron of a nucleus, should only be applied within a model and 
should not be used to compare among different models. Indeed, the Z271v model predicts 
the weakest (strongest) binding energy for the last neutron orbital in ^°Ca (^^O) among all 
the models. 

A case of particular interest is the doubly magic nucleus ^^O (upper- right-hand panel of 
Fig. HI) which appears to be particle unstable in spite of the many theoretical predictions 



to the contrary j20|, |2l|, |22, |23| . While the case of ^^O in particular, and the whole isotopic 
chain in general, deserves special attention and thus a separate publication, suffices to say 
that in all the self-consistent relativistic mean-field models considered here, the IcP^"^ orbital 
in ^^O is predicted to be bound by at least 2 MeV. 

We finish the discussion of Fig. 0] by addressing the model dependence of the neutron drip 
lines in calcium and zirconium. Predictions for the least bound neutron in the neutron-rich 
nuclei ™Ca (lower-left-hand panel) and ^^^Zr (lower-right-hand panel) are displayed as a 
function of the neutron skin in ^"'^Pb. Note that a "binding energy" of MeV indicates 
that the neutron is unbound and drips. The NL3 model predicts the last occupied neutron 
orbital in both nuclei {Ig^^'^ and 3p^/^, respectively) to be bound, albeit only weakly. This is 
in contrast to the S271 model for which both single-particle orbitals are bound but only for 
those parameterizations with a large neutron skin in ^''^Pb; that is, with a stiff symmetry 
energy. The Z271 model predicts both nuclei to be particle unstable for all values of Rn—Rp 
in 208pb. 

We conclude this section with a brief comment on the impact of a 1% measurement of the 
neutron radius in ^°^Pb on the neutron radius of other heavy nuclei that have been identified 
as promising candidates for atomic parity-violation experiments: barium, dysprosium, and 
ytterbium j2J, 2^ 2^ . Part of the appeal of these atoms is the existence of very close (nearly 



degenerate) levels of opposite parity that considerably enhance parity-violating amplitudes. 
Unfortunately, parity-violating matrix elements are contaminated by uncertainties in both 
atomic and nuclear structure. A fruitful experimental strategy for removing the sensitivity 
to the atomic theory is to measure ratios of parity-violation observables along an isotopic 
chain. This leaves nuclear-structure uncertainties, in the form of the neutron radius, as 
the limiting factor in the search for physics beyond the standard model 27, |2^ 2^. All 



three elements, barium, dysprosium, and ytterbium, have long chains of naturally occurring 
isotopes. While the experimental strategy demands a precise knowledge of neutron radii 
along the complete isotopic chain, we only correlate here (as means of illustration) the 
neutron radius of ^osp^ to the neutron radius of a member of the isotopic chain having 
a closed neutron shell (or subshell); we only consider here: ^^^Ba(Z = 56; iV = 82) j30l |. 
i58Dy(Z = 66; = 92), and ^'^'^Yh{Z = 70;N= 106). Note that for the open proton shell a 
spherical average is performed. That is, the factor of 2jk + 1 in Eq. ()13|) is simply replaced 
by the actual number of protons in the orbital. 

The neutron skins of ^^^Ba, ^^^Dy, and ^^^Yb, are correlated to the corresponding neutron 
skin of 208p]^ the three panels of Fig. We observe a tight linear correlation that is 
largely model independent. The linear regression coefficients (slope m and intercept b) have 
been enclosed in parenthesis. A theoretical spread of approximately 0.3 fm in the neutron 



radius of ^"J^p]^ ^g_g estimated in Refs. 27, \3m on the basis of best-fit models. Most of 



this spread is driven by the difference between relativistic and nonrelativistic models, which 
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FIG. 5: Skin-skin correlations for three heavy nuclei of possible relevance to atomic parity violation 
(^^^Ba, ^^'^Dy, and ^^^Yb) as a function of the neutron skin of ^'^^Pb for the three models considered 
in the text. Quantities in parenthesis represent linear regression coefficients (slope and intercept). 



has recently been attributed to the poorly known density dependence of the symmetry 
energy jl^. With the culmination of the the Parity Radius Experiment at the Jefferson 
Laboratory , the theoretical spread will be replaced by a genuine experimental error that 
is five times smaller, that is, Ai?„(^*'^Pb) ^ 0.056 fm. This 1% measurement of the neutron 
radius in ^''^Pb translates into a neutron radius uncertainty of Ai?„(^^^Ba) ^ 0.045 fm, 
ARX^^'Dj) ~ 0.034 fm, and AR^'^^^^) - 0.052 fm, respectively. 
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IV. CONCLUSIONS 



The neutron radius of ^"'^Pb, an observable highly sensitive to the density dependence 
of the symmetry energy, is correlated to ground-state observables of neutron-rich nuclei. 
We find that models with small neutron skins in ^osp]^ predict weak binding energies for 
the valence neutron orbitals in neutron-rich nuclei. Thus, models with the softest symmetry 
energy are the first ones to drip neutrons. Further, a tight correlation was found between the 
neutron skin of ^°^Pb and the neutron radius of a variety of elements (barium, dysprosium, 
and ytterbium) of possible relevance to atomic parity violation. 

The softening of the symmetry energy, which generates a thinner neutron skin in ^°^Pb, 
is accomplished by introducing additional (isoscalar-isovector) terms into the effective La- 
grangian. The new terms represent an important addition to the relativistic mean-field 
models, as they allow modifications to the neutron skin of ^osp]^ without compromising 
their success in reproducing a variety of ground-state observables. At first, it might seem 
surprising that the neutron radius of ^°^Pb, one of the most studied nuclei both theoreti- 
cally and experimentally, should be so poorly known. Indeed, theoretical estimates place 
an uncertainty in the neutron radius of ^'^'^Pb at about 0.3 fm. This is in contrast to its 
charge radius which is known — experimentally — to exquisite accuracy. The reason for this 
mismatch is twofold. First, electron scattering experiments, arguably the cleanest probe 
of nuclear structure, are only sensitive to the proton distribution. Second, best-fit models 
constrained to reproduce a variety of nuclear observables, such as charge densities, binding 
energies, and single-particle spectra, still predict a wide range of neutron radii for ^^'^Pb. 

Fortunately, the measurement of the neutron radius of ^ospb seems within reach. Indeed, 
the Parity Radius Experiment at the Jefferson Laboratory aims to measure the neutron 
radius in ^osp]^ accurately (at the 1% level) and model independently via parity-violating 
electron scattering. Such a measurement seems vital, as knowledge of a single isovector 
observable is sufficient to place stringent constraints on the model dependence of symmetry 
energy. Further, such a measurement will shed light on a variety of rich nuclear phenomena, 
ranging from the structure of neutron-rich nuclei to the structure of neutron stars. 

This work was supported in part by DOE grant DE-FG05-92ER40750. 
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